Fractional variational calculus in terms 
of a combined Caputo derivative 
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1. INTRODUCTION 

The history of Fractional Calculus (FC) goes back more 
than three centuries, when in 1695 the derivative of order 
a = 1/2 was described by Leibniz. Since then, the new 
theory turned out to be very attractive to mathematicians 
and many different forms of fractional operators were 
introduced: the Grunwald-Letnikov, Riemann-Liouville, 
Riesz, and Caputo fractional derivatives (see, e.g., Kilbas 
et al. [2006], Podlubny [1999], Samko et al. [1993]), and 
the more recent notions of Klimek [2002], Cresson [2007], 
and Jumarie [2007, 2010]. Besides mathematics, fractional 
^ derivatives and integrals appear in physics, mechanics, en- 
. gineering, elasticity, dynamics, control theory, economics, 
'biology, chemistry, etc. (cf. Baleanu [2009], Bastos et al. 
[2010], Mozyrska and Torres [2010], Rabei et al. [2007] and 
references therein) . The FC is nowadays covered by several 
books (e.g., Hilfer [2000], Kilbas et al. [2006], Miller and 
Ross [1993], Oldham and Spanier [1974], Podlubny [1999], 
Samko et al. [1993]) and a large number of relevant pa- 
pers (see, e.g., Agrawal [2007], Almeida and Torres [2010], 
Baleanu et al. [2010], Frederico and Torres [2008], Jumarie 
[2007, 2010], Klimek [2002], Rabei and Ababneh [2008], 
Tarasov [2008]). 

The calculus of variations is an old branch of optimization 
theory that has many applications both in physics and 
geometry (cf. Malinowska and Torres [2010b], Plakhov and 
Torres [2005], Troutman [1996], van Brunt [2004] and ref- 
erences therein). Apart from a few examples known since 
ancient times such as Queen Dido's problem (reported in 
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The Aeneid by Virgil), the problem of finding optimal 
curves and surfaces has been posed first by physicists such 
as Newton, Huygens, and Galileo. Their contemporary 
mathematicians, starting with the Bernoulli brothers and 
Leibniz, followed by Euler and Lagrange, invented the 
calculus of variations of a functional in order to solve those 
problems. Fractional Calculus of Variations (FCV) unifies 
the calculus of variations and the fractional calculus, by in- 
serting fractional derivatives into the variational integrals. 
This occurs naturally in many problems of physics or me- 
chanics, in order to provide more accurate models of phys- 
ical phenomena (see, e.g., Almeida et al. [2010], Atanack- 
ovic et al. [2008] , Baleanu [2009]). The FCV started in 1996 
with the work of Riewe [1996]. Riewe formulated the prob- 
lem of the calculus of variations with fractional derivatives 
and obtained the respective Euler-Lagrange equations, 
combining both conservative and nonconservative cases. 
Nowadays the FCV is a subject under strong research. 
Different definitions for fractional derivatives and integrals 
are used, depending on the purpose under study. Investi- 
gations cover problems depending on Riemann-Liouville 
fractional derivatives (see, e.g., Agrawal [2002], Almeida 
and Torres [2009], Atanackovic et al. [2008], Baleanu and 
Muslih [2005], El-Nabulsi and Torres [2008], Frederico and 
Torres [2007]), the Caputo fractional derivative (see, e.g., 
Agrawal [2007], Baleanu and Agrawal [2006], Malinowska 
and Torres [2010]), the symmetric fractional derivative 
(see, e.g., Klimek [2002]), the Jumarie fractional derivative 
(see, e.g., Almeida et al. [2010], Jumarie [2007, 2010]), 
and others (see, e.g., Agrawal [2007b], Cresson [2007], El- 
Nabulsi and Torres [2007], Frederico and Torres [2010]). 
Although the literature of FCV is already vast, much 
remains to be done. 



In this paper we extend the notion of the Caputo fractional 
derivative to the fractional derivative "-^D^'P, which is a 
convex combination of the left Caputo fractional derivative 
of order a and the right Caputo fractional derivative of 
order /3. This idea goes back at least as far as Klimek 
[2002], where based on the Riemann-Liouville fractional 
derivatives the symmetric fractional derivative was intro- 
duced. Klimek's approach is obtained choosing our pa- 
rameter 7 to be 1/2. Although the symmetric fractional 
derivative of Riemann-Liouville introduced by Klimek is 
a useful tool in the description of some nonconservative 
models, this type of differentiation does not seems suit- 
able for all kinds of variational problems. The hypothesis 
that admissible trajectories y have continuous symmetric 
fractional derivatives implies that y{a) = y{b) = (cf. 
Ross et al. [1994/95]). Therefore, the advantage of the 
fractional derivative D"'^ via Caputo lies in the fact that 
using this derivative we can describe a more general class 
of variational problems. It is also worth pointing out that 
our fractional derivative '~' D"'^ allows to generalize the 
results presented in Agrawal [2007]. 

The text is organized as follows. Section 2 presents some 
preliminaries. In Section 3 we introduce the fractional 
derivative D'^'P and provide the necessary concepts and 
results needed in the sequel. Our main results are stated 
and proved in Section 4. The fractional variational prob- 
lems under our consideration are formulated in terms 
of the fractional derivative "-^ D"'^ . We discuss the fun- 
damental concepts of a variational calculus such as the 
Euler-Lagrange equations for the basic (Subsection 4.1) 
and isoperimetric (Subsection 4.3) problems, as well as 
transversality conditions (Subsection 4.2). We end with 
conclusions in Section 5. 



2. PRELIMINARIES 

In this section we review the necessary definitions and facts 
from fractional calculus. For more on the subject we refer 
the reader to Kilbas et al. [2006], Oldham and Spanier 
[1974], Podlubny [1999], Samko et al. [1993]. 

Let / G Li{[a,b]) and < a < 1. We begin with the 
left and the right Riemann-Liouville Fractional Integrals 
(RLFI) of order a of a function /. The left RLFI is defined 
as 



al^fix) 



1 



r(a) 

while the right RLFI is given by 
1 



{x-tr-^f{t)dt, xe[a,b], (1) 



.ib.n=^) 



r(a) 



it-x)"-\f{t)dt, xe[a,b], (2) 



where r(-) represents the Gamma function. Moreover, 
al^f = xlt f = ,/ if / is a continuous function. 

Let / G AC {[a, b]), where AC {[a, b]) represents the space of 
absolutely continuous functions on [a, 6] . Using equations 
(1) and (2), we define the left and the right Riemann- 
Liouville and Caputo derivatives as follows. The left 
Riemann-Liouville Fractional Derivative (RLFD) is given 

by 



aD^fix) 



1 



(a; - t)-"/(f)dt 



r(l — a) dx jjj 
^—all^"f{x), xG[a,6], 



(3) 



dx 



the right RLFD by 



r(l — a) dx J, J. 



{t - x)-°'f{t)dt 



-^]-^b.fi^)^ xG[a,6], 



(4) 



the left Caputo Fractional Derivative (CFD) is defined by 

'iDZfix) 



ruh)iy^-'^-"jrf^'^'' 



(5) 



al^ "^/(a;), xe [a,b]. 



and the right CFD by 



__i_^„-.)-.i^/«),, 



x-^b 



(6) 



— )/(x), xG[a,6], 



where a is the order of the derivative. 

The operators (l)-(6) are obviously linear. We now present 
the rule of fractional integration by parts for RLFI (see for 
instance Brunetti et al. [1993]). Let 0<Q:<l,p>l, 
g > 1, and 1/p + 1/q < 1 + a. If ,g G Lp{[a,b]) and 
f eLg{[a,b]), then 

rb i-b 

g{x)J^f{x)dx = / f{x),I^g{x)dx. (7) 

J a 

In the discussion to follow, we will also need the following 
formulae for fractional integrations by parts: 

rb 

g{x) '^D^f{x)dx 



I x—b 



= /Wx/r"5(x)|^^„+ / fix),D^gix)dx, 

•^" (8) 

g{x) ^Dlf{x)dx 



= ^f{x)al',-'^g{x)\lJ^+ / f{x)aD^g{x)dx. 

They can be easily derived using equations (3)-(6), the 
identity (7), and performing integration by parts. 

3. THE FRACTIONAL ^D^^/^ OPERATOR 

Let a, /3 G (0, 1) and 7 G [0, 1]. Motivated by the diamond- 
alpha operator used in time scales (see, e.g., Malinowska 
and Torres [2009], Sidi Ammi et al. [2008]), we define the 
fractional derivative operator '^' D"'^ by 

^'D-^^ = 7^^: + (1 - ifxD', , (9) 

which acts on / G AC {[a, b]) in the following way: 



^'D^^^fix) = 7^^:/W + (1 - if^o'lfix). 



Note that 
c r,°',li 



Cn^ 



Cria.p _ 



D^''^fix) ^ 'iD'ifix), '^D'^.'^fix) = ),Dlf{x) 



For f = [/i,...,/Ar] : [a,b] ^ M^, (TV G N) and /, G 
AC {[a, b]), i — 1, . . . ,N, we have 

The operator (9) is obviously hnear. Using equations (8) it 
is easy to derive the following rule of fractional integration 
by parts for ^D^'f^: 

''" g{x) ^'D'^-^fix)dx ^ 7 [/(x)./i-"5(x)] Zl 



fix)Dt%gix)dx, (10) 



where -Dfl" = (1 — ^)aD^ + j xD^, which acts on / as 
DtVi^) = (1 - l)aD^.f{x)+j.D^f{x) . 

Let D denote the set of all functions y : [a, b] -^ R^ such 
that D"'' y exists and is continuous on the interval [a, b]. 
We endow D with the following norm: 



|y||i,. 



max |ly(x)|| + max ||'^i:'"''^y(x) 



i<x<h a<x<b 

where || • || stands for a norm in R^. 

Along the work we denote by diK, i = 1, . . . , M (M G N), 
the partial derivative of function K : R*^ -^ M with respect 
to its ith argument. 

Let A G M''. For simplicity of notation we introduce the 
operators [y] and {y}\ defined by 

[y]{x) = {xM^),''D^M^)) , 
{y}A(x) = {x,y{x),^D^^^y{x),Xu...,Xr) • 

4. CALCULUS OF VARIATIONS VIA ^D^'l^ 

We are concerned with the problem of finding minima (or 
maxima) of a functional ^7 : 2? — t- R, where P is a subset of 
D. The formulation of a problem of calculus of variations 
requires two steps: the specification of a performance 
criterion and the statement of physical constraints that 
should be satisfied. 

A performance criterion JT', also called cost functional (or 
cost), must be specified for evaluating the performance of 
a system quantitatively. We consider the following cost: 

J(y)= f L[y]{x)dx, 

J a 

where x G [a, b] is the independent variable, often called 
time; y(x) G R^ is a real vector variable, the functions 
y are generally called trajectories or curves; '^D"'^y(x) G 
R^ stands for the fractional derivative of y(a;); and L G 
C^{[a, b] X M^A^. R) is called a Lagrangian. 

Enforcing constraints in the optimization problem reduces 
the set of candidate functions and leads to the following 
definition. 

Definition 1. A trajectory y G D is said to be an ad- 
missible trajectory provided that it satisfies all of the 
constraints along interval [a,b]. The set of admissible tra- 
jectories is defined as 

2? := {y G D : y is admissible}. 



Wc now define what is meant by a minimum of J on T). 

Definition 2. A trajectory y G P is said to be a local 
minimizer for JT" on 2? if there exists (5 > such that 
J{y) < J{y) for ah y G 2? with ||y - y||i,oo < S. 

The concept of variation of a functional is central to the 
solution of problems of the calculus of variations. 

Definition 3. The first variation of JZ at y G D in the 
direction h G D is defined as 

A -7^ 1,^ r :^(y + eh) - J{y) 8 

6J{y; h := lim = -^ J(y + eh) 

provided the limit exists. 

Definition 4- A direction h G D, h 7^ 0, is said to be an 
admissible variation at y G 2? for J if 

(i) SJ'{y; h) exists; and 
(ii) y + eh G 2? for all sufficiently small e. 

The following well known result (see, e.g., [Troutman, 
1996, Proposition 5.5]) offers a necessary optimality con- 
dition for problems of calculus of variations based on the 
concept of variations. 

Theorem 5. (Troutman [1996]). Let ^J be a functional de- 
fined on 2?. Suppose that y is a local minimizer for Jj on 
v. Then, 5J{y; h) = for each admissible variation h in 

y- 

4-1 Elementary problem of the '~^ D'^P calculus of variations 



Let us begin with the following problem: 

J{y) = / L[y]{x)dx — > min 

J a 

over all y G D satisfying the boundary conditions 



y(a)=y^ y(6) - y^ y^y'' G 



cW 



(11) 



(12) 



Next theorem gives the fractional Eulcr -Lagrange equa- 
tion for the problem (11)-(12). 

Theorem 6. Let y = {yi, . . . ,yjq) be a local minimizer 
to problem (11)~(12). Then, y satisfies the system of N 
fractional Euler-Lagrange equations 



^,L[y]{x)+D^^'^^^N+^L[y]{x) 
for all x G [a, h\. 



0, 



.TV, (13) 



Proof. Suppose that y is a local minimizer for J . Let 
h be an arbitrary admissible variation for problem (11)- 
(12), i.e., hi{a) = hi{b) = 0, i = l,...,iV. Based on 
the differentiability properties of L and Theorem 5, a 
necessary condition for y to be a local minimizer is given 

by 

|:^(y + eh) =0, 



de 



e=0 



that is 

rb 



N+l 

^ d^L[y]{x)h^^i{x) 
.1=2 

N+l 

+ J2 ^N+^L[y]{xfD^^^h,^l{x) 



dx = 0. (14) 



Using formulae (10) for integration by parts in the second 
term of the integrand function, we get 



N+1 



i=2 



+ 7 



Y,^^L[y]{x)+Dfi^^^^M+^L[y]{x) 

N+l 

i=2 
'N+l 



hi-i{x)dx 

=b 



(1-7) 



x—a 
x—b 



= 0. 



Since hi(a) = hi{b) =0, i = 1, . . . ,N, by the fundamental 
lemma of the calculus of variations we deduce that 

a,i[y](x) + Dt"^^N+^L[y]{x) =0, i = 2,...,N + l, 
for all X e [a, b]. 

Observe that if a and /3 go to 1, then ^-D^. and a^*" can 

be replaced with ^; and ^Z?^ and ^Db with — ^ (sec, 
e.g., Podlubny [1999]). Thus, if 7 = 1 or 7 = 0, then for 
a, /3 — > 1 we obtain a corresponding result in the classical 
context of the calculus of variations (see, e.g., Troutman 
[1996]). 



4-2 ^ D!^'P transversality conditions 

Let I € {1,...,N}. Assume that y(a) = y°, yi{b) = y\, 
i = 1,...,A^ , i 7^ Z, but yi{b) is free. Then, hi{b) is free 
and by equations (13) and (15) we obtain 



lJl-^di+,L[y]{x) 



^{l-^)Jl-''dM+i+iL[y]{x) 



0. (16) 



Let us consider now the case when y(a) — y", yi{b) — y'l, 
i = 1,...,N , i 7^ I, and yi{b) is free but restricted by 
a terminal condition yi{b) < y'j^. Then, in the optimal 
solution y we have two possible types of outcome: yi{b) < 
yf or yi{b) — yf. li y{b) < yi, then there are admissible 
neighboring paths with terminal value both above and 
below yi{b), so that hi{b) can take either sign. Therefore, 
the transversality conditions is 

-fJl-''di+iL[y]{x) 

-(l-7)a/^^a„+l+^L[y](x)l =0 (17) 

- x—b 

for yi{b) < y^. The other outcome yi{b) ~ y^ only admits 
the neighboring paths with terminal value y{b)i < yi{b). 
Assuming, without loss of generality, that hi{b) > 0, this 
means that e < 0. Hence, the transversality condition, 
which has it root in the first order condition (14), must 
be changed to the inequality. For a minimization problem, 
the < type of inequality is called for, and we obtain 



jxll-''di+iL[y]{x) 

-(l-7)a/^^9A^^ 



HL[y]{x) 



X—b 



< (18) 



for yi{b) = y\ . Combining (17) and (18), we may write 
the following transversality condition for a minimization 
problem: 



^jl-''di+^L[y]{x) 



-{l--l)Jl-"dN+i+iL[y]{x) <0, 

- x^b 

vi{b)<yl 

{yi{b)-y\)[iJl-"di+,L[y]{x) 
-(l-7)a/^^a^^+l+^L[y](x) 



= 0. 



v^'^-' 4-3 The D"'^ isoperimetric problem 



Let us consider now the isoperimetric problem that con- 
sists of minimizing (11) over all y G D satisfying r isoperi- 
metric constraints 

cb 



Q'{y) 



G^[y]{x)dx^l,, j-1,. 



(19) 



, r, and boundary 



where G^' e Ci([a, fe] xR2^;M), j = 1, 
conditions (12). 

Necessary optimality conditions for isoperimetric problems 
can be obtained by the following general theorem (see, e.g., 
[Troutman, 1996, Theorem 5.16]). 

Theorem 7. Let J ,Q^ ,. . . ,Q^ be functionals defined in a 
neighborhood of y and having continuous first variations 
in this neighborhood. Suppose that y is a local minimum 
of (11) subject to the boundary conditions (12) and the 
isoperimetric constrains (19). Then, either: 
(i) for ah h^ G D, j = 1, . . . ,r, 

5g\yy) 5g\y-\i^) ■■■ 5g\y-W) 

sg^yX) '5e'(y;h2) •■• <5e'(y;h'^) 



5g^{y-\v') W{y.\y^ 



6g-{y-\f 



= 



(20) 



(ii) there exist constants Aj G M, j = 1, . . . , r, for which 

m 

5J{y: h) = ^ A,<5g-'(y; h) Vh e D. 

Note that condition (ii) of Theorem 7 can be written in 
the form 



5| J(y;h)-^A,g^(y;h) 



(21) 



Suppose now that assumptions of Theorem 7 hold but 
condition (i) docs not hold. Then, equation (21) is fulfilled 
for every h G D. Let us consider function h such that 
h(a) = h(&) = and denote by T the functional 



3 = 1 



Then, we have 
= 5T(y; h) 



de 



Hy + ^h) 



£ = 



'N+l 



J29^F{y}xix)h,^l{a 



N+l 



Y,^N+^F{y}^{xfD^'^h,-,i3 



dx, 



where the function F 



[a, b] X M2W X : 



is defined by conclusion that there exist constants Xj G M, j = 1, . 



F{y},{x) = L[y]{x) -J2x,Gny]{x). 

On account of the above, and similarly in spirit to the 
proof of Theorem 6, we obtain 



d,F{y}x{x) + Df■:^^^N+^F{y}x{x) =0, * = 2, ... TV. 

(22) 

Therefore, we have the following necessary optiniality 
condition for the isoperinietric problem: 

Theorem 8. Let assumptions of Theorem 7 hold. If y is a 
local minimizcr to the isoperimetric problem (11), (12) and 
(19), and condition (20) does not hold, then y satisfies the 
system of N fractional Euler-Lagrange equations (22) for 
all X e [a,b]. 

Suppose now that constraints (19) are characterized by 
inequalities 

G'iy)^ I G^[y]{x)dx<l,, j = i,...,r. 



dx+ {4>j{x)fdx = 0, 



In this case we can set 

'' 'G^[y]{x)--^ 
— a ^ 

J = 1, . . . , r, where (pj have the some continuity properties 
as Hi. Therefore, we obtain the following problem: 

J{y)^ / L{x,y{x),'' D'^'^y{x),cj,{x))dx ^ xmn 

J a 

where 4>{x) — [(/)i(x), . . . , (prix)], subject to r isoperimetric 
constraints 



G'[y]{x)- 



U 



{cP,{x)f 



dx — 0, j — 1, . 



and boundary conditions (12). Assuming that assumptions 
of Theorem 8 are satisfied, we conclude that there exist 



constants A, G 



J 



1, . . . , r, for which the system of 



equations 

Dt"^^N+^F{x, y(x), ""D^'^yix), Ai(x), . . . , A,(x), 0(x)) 

+ d,F{x, y(x), ^D^'^y{x), X,{x), . . . , A,(x), 0(x)) = 0, 

(23) 

i = 2,...,iV + l, where F = L + ^^=1 ^j (G' - ^ + <(']) 



and 



Xj(j>j{x)^o, i = 1, 



(24) 



hold for all x G [a, 6]. Note that it is enough to assume that 
the regularity condition holds for the constraints which are 
active at the local minimizer y. Indeed, suppose that I < r 
constrains, say Q^ , . . . ,G^ for simplicity, are active at the 
local minimizer y, and 

sg'irX) sg\y;h') ■■■ sg^yM) 

Sg'(y;h') Sg'(y;h') ■■■ Sg'{yM) 



sg'iyX) sgHy-y 



sg'iy-M 



^0 



for (independent) h^ G D, j = 1, . . . , L Since the inequality 
constraints tj'+i, ■ ■ ■ ,g^ are inactive, the conditions (24) 
are trivially satisfied by taking Xi+i ^ . . . — Xr ~ 0. On 
the other hand, since the inequality constraints g^ , . . . ,g^ 
are active and satisfy a regularity condition at y, the 



such that (23) holds follow from Theorem 8. Moreover, 
(24) is trivially satisfied for j = 1, . . . , I. 

5. CONCLUSION 

The FC is a mathematical area of a currently strong 
research, with numerous applications in physics and engi- 
neering. The fractional operators are non-local, therefore 
they are suitable for constructing models possessing mem- 
ory effect. This gives several possible applications of the 
FCV in describing non-local properties of physical systems 
in mechanics or electrodynamics. In this note we extend 
the notions of Caputo fractional derivative to the fractional 
derivative ^ D"'^ . We emphasize that this derivative allows 
to describe a more general class of variational problems 
and, as a particular case, we get the results of Agrawal 
[2007]. 

Knowing the importance and relevance of multiobjective 
problems of the calculus of variations in physics and 
engineering, our further research will continue towards 
multiobjective FCV. This is a completely open research 
area and will be addressed elsewhere. 
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